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Summary. — An overview is presented of the fundamentals of continuum-polaron 
physics, which provide the basis of the analysis of polaron effects in ionic crystals 
and polar semiconductors. The present paper deals with "large", or "continuum", 
polarons, as described by the Frohlich Hamiltonian. The emphasis is on the polaron 
optical absorption. 
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1. The Polaron Concept 

As is generally known, the polaron concept was introduced by Landau in 1933 1 . 
Gradually a distinction was made between polarons in the continuum approximation 
where long-range electron-lattice interaction prevails ("Frohlich" polarons) and polarons 
for which the short-range interaction is essential (Holstein, Holstein- Hubbard models). 
After a period of initial theoretical 012131011311113111 an d experimental |5] works, a 
first international meeting was devoted to the subject in St. Andrews (1965, 10 ). Ten 
years later a second international meeting in which the progress in the field of polarons 
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realized by 1972 was recorded - took place in Antwerpen Among the review papers 
and books covering the subject, we mention Refs. [T^l ITTHfT^l fTKj . 

Although the present course in Varenna 2005 mainly deals with "small" polarons (Hol- 
stein, Holstein-Hubbard) , my contribution is devoted to "large" polarons. The reader is 
referred to most other papers in the present book for an overview of the properties of 
Holstein related polarons. In the present volume the contributions ^3 El El d!| a l so 
deal with the Frohlich polaron. 

Let us recall the polaron concept for newcomers. A charge placed in a polarizable 
medium is screened. Dielectric theory describes the phenomenon by the induction of a 
polarization around the charge carrier. The induced polarization will follow the charge 
carrier when it is moving through the medium. The carrier together with the induced 
polarization is considered as one entity (see Fig.^). It was Pekar who came up with 
the term polaron |21j . The physical properties of a polaron differ from those of a band- 
carrier. A polaron is characterized by its binding (or self-) energy Eq, an effective mass 
to* and by its characteristic response to external electric and magnetic fields (e. g. dc 
mobility and optical absorption coefficient). 
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Fig. 1. - Artist view of a polaron. A conduction electron in an ionic crystal or a polar semicon- 
ductor repels the negative ions and attracts the positive ions. A self-induced potential arises, 
which acts back on the electron and modifies its physical properties. (From |2U|.) 

If the spatial extension of a polaron is large compared to the lattice parameter of the 
solid, the latter can be treated as a polarizable continuum. This is the case of a large, or 
continuum (Frohlich) polaron. When the self-induced polarization caused by an electron 
or hole becomes of the order of the lattice parameter, a small (Holstein) polaron can arise. 
As distinct from large polarons, small polarons are governed by short-range interactions. 

IT. The Frohlich Hamiltonian. - Frohlich proposed a model Hamiltonian for the 
"large" polaron through which its dynamics is treated quantum mechanically ( "Frohlich 
Hamiltonian"). The polarization, carried by the longitudinal optical (LO) phonons, is 
represented by a set of quantum oscillators with frequency wlo, the long- wavelength 
LO-phonon frequency, and the interaction between the charge and the polarization field 
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is linear in the field [H]: 

(1) H =Y~+T, ^Lo4«k + ^(^ fe a k e lkr + ^*4e- ik ' r ), 

TOfc k k 

where r is the position coordinate operator of the electron with band mass nib, P is its 
canonically conjugate momentum operator; a k and ak are the creation (and annihilation) 
operators for longitudinal optical phonons of wave vector k and energy hu^o- The Vk 
are the Fourier components characterising the electron-phonon interaction 



(2) 
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The strength of the electron-phonon interaction is expressed by a dimensionless coupling 
constant a. which is defined as: 



(3) mb ' 1 ' 



In this definition, Soo and Eq are, respectively, the electronic and the static dielectric 
constant of the ionic crystal or polar semiconductor. 

In Ref. 1201) the Frohlich coupling constant is given for a few solids( 1 ). 

In deriving the form of 14, expressions © and ©, it was assumed that (i) the 
spatial extension of the polaron is large compared to the lattice parameter of the solid 
("continuum" approximation), (ii) spin and relativistic effects can be neglected, (iii) the 
band-electron has parabolic dispersion, (iv) in line with the first approximation it is also 
assumed that the LO-phonons of interest for the interaction, are the long-wavelength 
phonons with constant frequency wlo 

The model, represented by the Hamiltonian l[T]l. (which up to now could not been 
solved exactly) has been the subject of extensive investigations, see, e. g., Refs. [31 ITT71 
1121 II 11 1221 1131 ITU I23| . In what follows the key approaches of the Frbhlich-polaron theory 
are briefly summarized with indication also of their relevance for the polaron problems 
in nanostructures. 



2. The Structure of the Polaron Excitation Spectrum and Optical Absorp- 
tion by Polarons 

2T. Optical absorption at weak coupling. - At zero temperature and in the weak- 
coupling limit, the optical absorption is due to the elementary polaron scattering process, 
schematically shown in Fig.[2J 



( 1 ) In some cases, due to lack of reliable experimental data to determine the electron band 
mass, the values of a are not well established. 
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Fig. 2. - Elementary polaron scattering process describing the absorption of an incoming photon 
and the generation of an outgoing phonon. 



In the weak-coupling limit (a<l) the polaron absorption coefficient was first studied 
by V. Gurevich, I. Lang and Yu. Firsov [21]. Their optical-absorption coefficient is 
equivalent to a particular approximation of the result of J. Tempere and J. T. Devreese 
(Ref. [23 ) j with the dynamic structure factor S(k, il) corresponding to the Hartree-Fock 
approximation (see also 120], p- 585). 

The result of Ref. [23] gives the Frohlich polaron optical absorption coefficient for a 
many-polaron gas (in terms of the dynamic structure factor S*(k, 51)) exactly to order a. 
The real part of the optical conductivity of the 7V-polaron system for any N is obtained 
in Ref. |25| in an intuitively transparent form 



where n — N/V is the density of charge carriers and k x is the ^-component of the 
wave vector (the electric field is along the x-axis). This approach to the many-polaron 
optical absorption allows one to include the many-body effects in terms of the dynamical 
structure factor 



of the electron ( or hole ) system. Substituting Eq. (01 in Eq. J3J , the following analytical 
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result was obtained in Ref. |25|f 2- ): 

+00 

(6) Re[<r(fi)] = ne 2 ^-=^— / dkk 2 S(k,n - w LO ). 

O \f2lT\l i J 


For N = 1, substituting a (5- function for the dynamical structure factor 

ftfc 2 



S(k,ft) 



2™^ 



- (ft - UJ LO ) 



in the general expression J2J, the one-polaron limit Ref. is retrieved. This one- 

polaron result is derived 23 by considering a process, in which the initial state consists 
of a photon of energy htt and a polaron in its ground state, and the final state consists 
of an emitted LO phonon with energy fko^o and the polaron, scattered into a state with 
momentum k and kinetic energy (hk) 2 / (2mb) = h(fl — cjlo)- The many-polaron result 
Q is a generalization of the one-polaron picture. The contribution which corresponds 
to the scattering of a polaron into the momentum state k and energy h(£l — lulo) is now 
weighed by the dynamical structure factor S(k, ft — cjlo) of the electron (or hole) gas. 

In the zero-tempcrature limit, within the weak coupling approximation, the real part 
of the polaron optical conductivity has the following analytic expression (see |27|) : 

o „ -> 3/2 

' . a 



(7) Re a (ft) = -—^(5 (O) + Q l° 3 V" ~ ^ LO e (Q - lu lo ) , 

2m* 6ml, il a 



where 



e(ft-WLo) 



1 if ft > u LO , 
if ft < wlo 



The spectrum of the real part of the polaron optical conductivity Q is represented in 
Fig. |3 The absorption coefficient for absorption of light with frequency ft by a gas of 
free polarons for a — > takes the form 



(8) r f(fl )=) = Jli(fl) + i^JI^ e ( ft- WLO ), 

e cn 2m* e cn 3mbil J V ^lo 

where eo is the dielectric permittivity of the vacuum, n is the refractive index of the 
medium, no is the concentration of polarons. A simple derivation in Ref. |28j using a 
canonical transformation method gives the absorption coefficient of free polarons, which 
coincides with the second term in the result JSJ- A step function in (JSJ reflects the fact 
that at zero temperature the absorption of light accompanied by the emission of a phonon 



( 2 ) A denominator of the prefactor in the Eq. (23) of Ref. |25| contain a misprint: there should 
be \p2 instead of 2. 
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can occur (i) at £1 = and (ii) if the energy of the incident photon is larger than that 
of a phonon (Cl > cJlo)- The second term in the absorption coefficient JHJ was obtained 
in Ref. |24| in the low-concentration region. In the weak-coupling limit, according to 
Eq. |JSJ, the absorption spectrum consists of a central peak [oc <5(f2)] and a "one-phonon 
sideband" . At nonzero temperature, the absorption of a photon can be accompanied not 
only by emission, but also by absorption of one or more phonons. 
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Fig. 3. - Polaron optical conductivity for a = 1 in the weak-coupling approximation. (After 
|27|. p. 92.) A S- like central peak (at Q = 0) is schematically shown by a vertical line. 

Experimentally, this one-phonon-sideband structure has been observed for free po- 
larons e.g. in the infrared absorption spectra of CdO-films, see Fig.0J In CdO, which 
is a weakly polar material with a ss 0.74, the polaron absorption band is observed in 
the spectral region between 6 and 20 /mi (above the LO phonon frequency). The differ- 
ence between theory and experiment in the wavelength region where polaron absorption 
dominates the spectrum is due to many-polaron effects. 

2'2. Optical Absorption at Strong Coupling. - The problem of the structure of the 
large polaron excitation spectrum constituted a central question in the early stages of 
the development of polaron theory. The exactly solvable polaron model of Ref. [3T] (see 
Appendix was used to demonstrate the existence of the so-called "relaxed excited 
states" of large polarons [32]. In Ref. [33], and after earlier intuitive analysis, this problem 
was studied using the classical equations of motion and Poisson-brackets. The insight 
gained as a result of those investigations concerning the structure of the excited polaron 
states, was subsequently used to develop a theory of the optical absorption spectra of 
polarons. The first work was limited to the strong coupling limit [33]. Ref. [33] is the 
first work that reveals the impact of the internal degrees of freedom of polarons on their 
optical properties. 

The optical absorption of light by free Frohlich polarons was treated in Ref. [33] using 
the polaron states obtained within the adiabatic strong-coupling approximation. It was 
argued in Ref. |33], that for sufficiently large a (a > 3), the (first) relaxed excited state 
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Fig. 4. - Optical absorption spectrum of a CdO-film with the carrier concentration n — 5.9 x 
10 cm" " 3 at T = 300 K. The experimental data (solid dots) of Ref. Hj] are compared to different 
theoretical results: with (solid curve) and without (dashed line) the one-polaron contribution 
of Ref. 12 II and for many polarons (dash-dotted curve) of Ref. |25| . The material parameters 
of CdO used for calculations: a — 0.74 1291 . ojlo = 490 cm -1 (from the experimental optical 
absorption spectrum, Fig. 2 of Ref. EH)), m b = 0.11m e |3I], e (0) = 21.9, e (oo) = 5.3 (HJ. 



(RES) of a polaron is a relatively stable state, which can participate in optical absorption 
transitions. This idea was necessary to understand the polaron optical absorption spec- 
trum in the strong-coupling regime. The following scenario of a transition, which leads 
to a "zero-phonon" peak in the absorption by a strong-coupling polaron, can then be 
suggested. If the frequency of the incoming photon is equal to Ores = 0.065q; 2 wlo, the 
electron jumps from the ground state (which, at large coupling, is well-characterized by 
"s" -symmetry for the electron) to an excited state ("2p"), while the lattice polarization 
in the final state is adapted to the "2p" electronic state of the polaron. In Ref. 
considering the decay of the RES with emission of one real phonon it is argued, that the 
"zero-phonon" peak can be described using the Wigner-Weisskopf formula valid when 
the linewidth of that peak is much smaller than hui^o ■ 

For photon energies larger than Ores + wlo, a transition of the polaron towards 
the first scattering state, belonging to the RES, becomes possible. The final state of 
the optical absorption process then consists of a polaron in its lowest RES plus a free 
phonon. A "one-phonon sideband" then appears in the polaron absorption spectrum. 
This process is called one-phonon sideband absorption. The one-, two-, ... K-, ... phonon 
sidebands of the zero-phonon peak give rise to a broad structure in the absorption spec- 
trum. It turns out that the first moment of the phonon sidebands corresponds to the 
Franck-Condon (FC) frequency f2pc = 0.141a 2 ct;LO- To summarize, the polaron optical 
absorption spectrum at strong coupling is characterized by the following features (at 
T = 0): 
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a) An intense absorption peak ( "zero-phonon line") appears, which corresponds to a 
transition from the ground state to the first RES at Ores- 

b) For f2 > Ores +wlo, a phonon sideband structure arises. This sideband structure 
peaks around fipc- 

The qualitative behaviour predicted in Ref. |34| . namely, an intense zero-phonon (RES) 
line with a broader sideband at the high-frequency side, was confirmed after an all- 
coupling expression for the polaron optical absorption coefficient had been studied |35| . 

2'3. Optical Absorption of Continuum Polarons at Arbitrary Coupling. - The polaron 
absorption coefficient r(fi) of light with frequency Vt at arbitrary coupling was first 
derived by the present author and collaborators (see also [US])- It was represented 
in the form 

(9) rm) = -J-^. 

ne c m b [Q _ Re £(ft)] 2 + [1m S(O)] 2 

This general expression was the starting point for a derivation of the theoretical optical 
absorption spectrum of a single Frohlich polaron at all electron-phonon coupling strengths 
in Ref. [35] • S(O) is the so-called "memory function", which contains the dynamics of the 
polaron and depends on f2, a and temperature. The key contribution of the work in |35) 
was to introduce r(f2) in the form and to calculate ReS(fi), which is essentially a 
(technically not trivial) Kramers-Kronig transform of the more simple function Im£(f2). 
The function lm£(0) had been formally derived for the Feynman polaron |37| who 
studied the polaron mobility fi from the impedance function, i. e. the static limit 



lim 



ImE(fi) 



fi-»o \ SI 

The nature of the polaron excitations was clearly revealed through this polaron optical 
absorption obtained in (351 l3l)| . It was demonstrated in [33] that the Franck-Condon 
states for Frohlich polarons are nothing else but a superposition of phonon sidebands. It 
was also established in |35| that a relatively large value of the electron-phonon coupling 
strength (a > 5.9) is needed to stabilise the relaxed excited state of the polaron. It was, 
further, revealed that at weaker coupling only "scattering states" of the polaron play a 
significant role in the optical absorption |351 \'6H\ . 

In the weak coupling limit, the optical absorption spectrum Q of the polaron is 
determined by the absorption of radiation energy, which is reemitted in the form of 
LO phonons. For a > 5.9, the polaron can undergo transitions toward a relatively 
stable RES (see Fig. [5J. The RES peak in the optical absorption spectrum also has 
a phonon sideband-structure, whose average transition frequency can be related to a 
FC-type transition. Furthermore, at zero temperature, the optical absorption spectrum 
of one polaron exhibits also a zero-frequency "central peak" [oc <5(fi)]. For non-zero 
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Fig. 5. - Optical absorption spectrum of a Frohlich polaron for a = 5 and a = 6. The RES 
peak is very intense compared with the FC peak. The frequency fl/tJLO = v is indicated by 
the dashed lines. The 5-like central peaks (at Q — 0) are schematically shown by vertical lines. 
(From Devreese et al. Ref. 1351 .1 



temperature, this "central peak" smears out and gives rise to an "anomalous" Drudc- 
type low-frequency component of the optical absorption spectrum. 

For example, in Fig. [S] (from Ref. 35 ), the main peak of the polaron optical ab- 
sorption for a = 5 at f2 = 3.51ojlo is interpreted as due to transitions to a RES. A 
"shoulder" at the low-frequency side of the main peak is attributed to one-phonon tran- 
sitions to polaron- "scattering states" . The broad structure centered at about f2 = 6.3wlo 
is interpreted as a FC band. As seen from Fig. [31 when increasing the electron-phonon 
coupling constant to a=6, the RES peak at fi = 4.3cjlo stabilizes. It is in Ref. (35] that 
the all-coupling optical absorption spectrum of a Frohlich polaron, together with the role 
of RES-states, FC-states and scattering states, was first presented. Based on Ref. [3*5] . 
it was argued that it is rather Holstein polarons that determine the optical properties 
of the charge carriers in oxides like SrTiOa, BaTiOa [35], while large weak coupling po- 
larons could be identified e. g. in CdO 0U]- The Frohlich coupling constants of polar 
semiconductors and ionic crystals are generally too small to allow for a static "RES" . 
In Ref. 01] the RES-peaks of Ref. j3U were involved to explain the optical absorption 
spectrum of Pr 2 Ni04.22- Further study of the spectra of is called for. The RES-like 
resonances in r(f2), Eq. ©, due to the zeros of — Re can effectively be displaced 

to smaller polaron coupling by applying an external magnetic field B, in which case the 
contribution for what is formally a "RES-type resonance becomes f2 — lo c — Re S(fi) = 
(u> c = eB/rribC is the cyclotron frequency). Resonances in the magnetoabsorption gov- 
erned by this contribution have been widely observed and analysed 021 to [SJ (see also 
Subsection 13 21 below) . 

Recent interesting numerical calculations of the optical conductivity for the Frohlich 
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Fig. 6. - Left-hand panel: Monte Carlo optical conductivity spectra of one Frohlich polaron 
for the weak-coupling regime (open circles) compared to the second-order perturbation theory 
(dotted lines) for a = 0.01 and a = 1 and to the analytical DSG calculations 1351 (solid lines). 
Right-hand panel: Monte Carlo optical conductivity spectra for the intermediate coupling regime 
(open circles) compared to the analytical DSG approach |35| (solid lines). Arrows point to the 
two- and three-phonon thresholds. (From Ref. |55|.) 



polaron performed within the Diagrammatic Quantum Monte Carlo (DQMC) method 
|55| . see Fig. El fully confirm the essential analytical results derived in Ref. for 
a < 3. In the intermediate coupling regime 3 < a < 6, the low-energy behavior and the 
position of the RES-peak in the optical conductivity spectrum of Ref. [55] follow closely 
the prediction of Ref. |35|. There are some minor qualitative differences between the 
two approaches in the intermediate coupling regime: in Ref. [S3], the dominant ("RES") 
peak is less intense in the Monte-Carlo numerical simulations and the second ("FC") 
peak develops less prominently. The following qualitative differences exist between the 
two approaches in the strong coupling regime: in Ref. |55| . the dominant peak broadens 
and the second peak does not develop, giving instead rise to a flat shoulder in the 
optical conductivity spectrum at a = 6. These qualitative differences from the optical 
absorption spectrum of Ref. |35| can be tentatively attributed to the optical processes 
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Table I. - Polaron parameters Mq, Mi, Mq obtained from the diagrammatic Monte Carlo results 



a 


M< MC) 


m »(MC) 


M< MC) 


M[ MC) /a 


P (MC) 
^0 


0.01 


0.00249 


1.0017 


1.5706 


0.634 


-0.010 


1 


0.24179 


1.1865 


1.5657 


0.65789 


-1.013 


3 


0.67743 


1.8467 


1.5280 


0.73123 


-3.18 


4.5 


0.97540 


2.8742 


1.5219 


0.862 


-4.97 


5.25 


1.0904 


3.8148 


1.5022 


0.90181 


-5.68 


6 


1.1994 


5.3708 


1.4919 


0.98248 


-6.79 


6.5 


1.30 


6.4989 


1.5417 


1.1356 


-7.44 


7 


1.3558 


9.7158 


1.5175 


1.2163 


-8.31 


8 


1.4195 


19.991 


1.4981 


1.3774 


-9.85 



with participation of two |56| or more phonons. The above differences can arise also due 
to the fact that, starting from the Feynman polaron model, one-phonon processes are 
assigned more oscillator strength and the RES of therefore tends to be more stable 
as compared to the Monte-Carlo result. The comparison of the analytical and numerical 
results for the optical absorption of the Frohlich polarons shown in Fig. [S]deserves further 
study. An independent numerical simulation is called for. 

2'4. Sum Rules for the Optical Conductivity Spectra. - In this section the sum rules for 
the optical conductivity spectra of Frohlich polarons obtained within the DSG approach 
[55] are compared with those found from the DQMC results [313- The values of the 
polaron effective mass for the Monte Carlo approach are taken from Ref. • In Tables 
[I] and m I show the polaron ground-state Eq and the following parameters calculated 
using the optical conductivity spectra: 



(10) M = J Rea(uj)duj, 

ujRecr(uj)duj, 

where u> max is the upper value of the frequency available from Ref. |55| . 
(12) M = — + / Recr(Lo)dLo, 

where m* is the polaron mass, the optical conductivity is calculated in units noe 2 /(mf,WLo); 
to* is measured in units of the band mass m^, and the frequency is measured in units of 
wlo- The values of u> max are: w max = 10 for a = 0.01, 1 and 3, u> max — 12 for a = 4.5, 
5.25 and 6, w max = 18 for a = 6.5, 7 and 8. 

The parameters corresponding to the Monte Carlo calculation are obtained using 
the numerical data kindly provided to the author by A. Mishchenko [57j. The optical 
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Table II. 


- Polaron parameters Mo, Mi, Mo 


obtained within the path-integral approach 


a 


M (DSG) 


^♦(Feynman) 


M (DSG) 


Af< DSG) /a 


i-i(Fcynman) 


0.01 

1 

3 

4.5 

5.25 

6 

6.5 

7 

8 


0.00248 

0.24318 

0.69696 

1.0162 

1.1504 

1.2608 

1.3657 

1.4278 

1.4741 


1.0017 
1.1957 

1.8912 
3.1202 
4.3969 
6.8367 
9.7449 
14.395 
31.569 


1.5706 
1.5569 
1.5275 
1.5196 
1.5077 
1.4906 
1.5269 
1.5369 
1.5239 


0.633 

0.65468 

0.71572 

0.83184 

0.88595 

0.95384 

1.1192 

1.2170 

1.4340 


-0.010 

-1.0130 

-3.1333 

-4.8394 

-5.7482 

-6.7108 

-7.3920 

-8.1127 

-9.6953 



conductivity obtained within the path-integral approach of Ref. [35) exactly satisfies the 
sum rule 



(13) 



2m* 



f°° IT 

I Re a (u>) dio = —. 
Ji 2 



Since the optical conductivity obtained from the Monte Carlo results [55] is known only 
within a limited interval of frequencies 1 < lu < w max , we calculate the integral in Eq. 
(|12fl for the DSG-approach over the same frequency interval as for the Monte Carlo 
results |55|. 

The comparison of the resulting zero frequency moments Mq MC ^ and Mq DSG ^ with 
each other and with the value ir/2 = 1.5707963... corresponding to the right-hand-side of 



M(MC)_^(DSG) 



on the interval a < 8 is 



the sum rule (|13fl shows that the difference 
smaller than the absolute value of the contribution of the "tail" of the optical conductivity 
for lo > w mal to the integral in the sum rule ifHUl : 



/•oo 

(14) / Ret7 (DSG) (u)du 



M 



>(DSG) 



Within the accuracy determined by the neglect of the "tail" of the part of the spectrum 
for uj > w mm , the contribution to the integral in the sum rule l|13|) f° r the optical 
conductivity obtained from the Monte Carlo results |55| agrees well with that for the 
optical conductivity found within the path-integral approach of Ref. |35j. Hence, the 
conclusion follows that the optical conductivity obtained from the Monte Carlo results 
J[55\j satisfies the sum rule within the aforementioned accuracy. 

We analyze also the fulfillment of the "LSD" polaron ground-state theorem introduced 
in [SH 



(15) 



E o (a)-E o {0) = / — / ojReo-(uj,a')duj 

7T Jo Ot' J 
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a 

Fig. 7. - Test of the ground-state theorem for a Frohlich polaron from Ref. |59| using different 
optical conductivity spectra, DSG from Ref. 1351 and MC from Ref. |55|. The notations are 
explained in the text. 



using the first-frequency moments m[ MC ^ and M i DSG ' ) . The results of this comparison 
are presented in Fig. [7| The dots indicate the polaron ground-state energy calculated 
using the Feynman variational principle. The solid curve is the value of Eq (a) calculated 
numerically using the optical conductivity spectra and the ground-state theorem with 
the DSG optical conductivity for a polaron, 

(16) 4 DSG) (^-- f — /" w Re CT ( DSG > (<<;,<*') 

Wo Jo 

(DSG) 

The dashed and the dot-dashed curves are the values obtained using M y ,> (a) and 
M[ MCS> (a), respectively: 

(17) E (a) = / —r I ujRea(uj,a')duj = / da \ ' . 

TT Jo Oi' J Q 7T J a' 

As seen from the figure, i?Q DSG ' 1 (a) to a high degree of accuracy coincides with the 
variational polaron ground-state energy. Both E^ DSG ^ (a) and _Eq MC - ) (a) differ from 

( DSC 1 1 

Eq (a) due to the integration over the same limited frequency interval. However, 
4 DSG) (a) and E^ C) (a) are very close to each other. Herefrom, the conclusion follows 
that for integrals over the finite frequency region characteristic for the polaron optical 
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absorption (i. e., omitting the "tails"), the function E^ 1 ^ (a) reproduces very well the 
function E^ SG ^ (a). 

From the comparison of £?q DSG ' (a) with E^ 1 ^ (a) it follows that the contribution to 
the integral in l|17[) with the limited frequency region, which approximates the integral in 
the right-hand side of the "LSD" ground state theorem l|15[l. for the optical conductivity 
obtained from the Monte Carlo results |55| agrees with a high accuracy with the corre- 
sponding contribution to the integral in (|17|l for the optical conductivity found within 
the path- integral approach of Ref. |35|. Because for the path- integral result, the inte- 
gral Sg DSG ' ) (a) noticeably differs from the available integral E^ 080 ' (a), a comparison 
between the Feynman polaron ground state energy Eq and the integral E^ SG ^ (a) is 
not representative. Similarly, a comparison between the polaron ground state energy 
obtained from the Monte Carlo results and the integral Eq MC ^ (a) would require to over- 
come the present limitation for the available optical conductivity spectrum |55j . 

In conclusion, the Monte Carlo optical conductivity spectrum for higher frequencies 
than w ma x of [55] is needed in order to check the fulfillment of the sum rules l|13[) and 
(|15|l with a higher accuracy. 

3. Scaling Relations for Frohlich Polarons 

Several scaling relations connect the polaron self-energy, the effective mass, the impedance 
Z and the polaron mobility /i in 2D to their counterpart in 3D. Those relations were intro- 
duced in Ref. jSHj (Peeters and Devreese), they are satisfied at the level of the Feynman 
polaron model and are listed here: 

, s 2 /3tt 
(18a) E 2D (a) = -E 3D ( —a 

(186) 



3 " V 4 



(18c) Z 2 u(a,Q) = Z 3V (^a,Q 

/37T 

(18(f) M2d(«) = M3D [ —£Ct 

In Eq. I|186(l . ("T-^rj) an d (w3d) are, respectively, the polaron and the electron- 
band masses in 2D (3D). Expressions i|18|) provide a straightforward link between po- 
laron characteristics in 3D with those in 2D. The scaling relations imply that polaron 
effects may be observable in 2 dimensions at substantially lower electron-phonon cou- 
pling than in 3 dimensions. Against this background, an extensive analysis of polaronic 
phenomena was undertaken for GaAs heterostructures, quantum wells, superlattices, 2D 
oxide-structures, shallow donors ("bound polarons") and also for quantum dots, see e. g. 
Refs. [H] to El and the reviews [T31I2U] . 
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Fig. 8. - Upper panel: Polaron free energy in 3D (squares) and 2D (open circles) obtained by 
TPC'2001 E3| for (3 = 10. The data for F SD (a) are interpolated using a polynomial fit to the 
available four points (dotted line) . Lower panel: Demonstration of the PD-scaling cf. Peeters- 
Devreese |6(J|: the polaron free energy in 2D obtained by TPC'2001 §55 ; for (3 — 10 (squares) is 
very close to the PD-scaled according to Peeters-Devreese |60| values of the polaron free energy 
in 3D from TPC'2001 for (3 = 10 (open triangles). The data for ^F SD (^) are interpolated 
using a polynomial fit to the available four points (solid line). 
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3'1. Check of the scaling relation for the path integral Monte Carlo results for the 
polaron free energy. - The fulfillment of the PD-scaling relation |6()| is now checked for 
the path integral Monte Carlo results for the polaron free energy. The path integral 
Monte Carlo results of Ref. for the polaron free energy in 3D and in 2D were given for 
a few values of temperature and for some selected values of a. For a check of the scaling 
relation, the values of the polaron free energy at (3 = 10 are taken from Ref. |^ in 3D 
(Table I of Ref. 65 , for four values of a) and in 2D (Table II of Ref. 65 , for two values 
of a) and plotted in Fig. El upper panel, with squares and open circles, correspondingly. 
In Fig. El lower panel, the available data for the free energy from Ref [HSj are plotted 
in the following form inspired by the l.h.s. and the r.h.s parts of Eq. $18a\) : F%d (a) 
(squares) and §1*3.0 (-^p)(open triangles). As follows from the figure, the path integral 
Monte Carlo results for the polaron free energy in 2D and 3D very closely follow the 
PD-scaling relation of the form given by Eq. \18a\) : 



3'2. Magnetoabsorption spectra of polarons. - The results on the polaron optical ab- 
sorption |35U27j paved the way for an all-coupling path-integral based theory of magneto- 
optical absorption by polarons (see Ref. 021 > F- Peeters and J. T. Devreese). This work 
was i. a. motivated by the insight that magnetic fields can stabilise the relaxed excited 
polaron states, so that information on the nature of relaxed excited states might be 
gained from the cyclotron resonance of polarons. A quantitative interpretation of the 
high-precision cyclotron resonance experiment on AgBr and AgCl in a range of magnetic 
fields up to 15 T 03] by the theory of Ref. 05| provided one of the most convincing 
and clearest demonstrations of polaron features in solids. For a review of the subsequent 
developments in the field of polaron cyclotron resonance, the reader is referred to Refs. 
@H1 EH! EH EH1 HIZI [US E3 an d references therein. The energy spectra of such polaronic 
systems as shallow donors ("bound polarons"), e. g., the D and D~ centres, constitute 
the most complete and detailed polaron spectroscopy realised in the literature, see for 
example Fig. El 

4. Many-Polaron Systems 

4T. Ground State and the Optical Absorption Spectra of Many-Polaron Systems. - 
The earliest studies on many polarons concerned their role in heavily doped polar semi- 
conductors HSU E3 IZD- The possibility that polarons and bipolarons play a role in 
high-T c superconductors has renewed interest in the physical properties of many-polaron 
systems and, in particular, in their optical properties. Theoretical treatments have been 
extended from one-polaron to many-polaron systems both of Frohlich polarons (see, for 
example, [721 EH CU EH1 E3 EH IZ3 CHI ) and of Holstein polarons (see, e. g., [23 EDI ED). 

For the weak-coupling regime, which is realized in most polar semiconductors, the 
ground-state energy of a gas of interacting large polarons was derived in Ref. [73] by 



(19) 
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Fig. 9. - The Is — ► 2p ± ,2p z transition energies as a function of a magnetic field for a donor 
in GaAs. We compare our theoretical results for the following cases: (a) without the effect 
of polaron and band nonparabolicity (thin dashed curves); (b) with polaron correction (dotted 
curves); (c) including the effects of polaron and band nonparabolicity (solid curves) to the 
experimental data of Ref. (solid dots). (From Ref. pb].) 



introducing a variational wave function: 

(20) \1>ldb) = U\<I>)\<p a ), 

where \<p e i) represents the ground-state many-body wave function for the electron (or 
hole) system, \(f>) is the phonon vacuum and U is a many-body unitary operator that 
determines a canonical transformation for a fermion gas interacting with a boson field: 

N 

(21) U = exp { E (/ka k e* k ^ - fZa+e~ ikr 

J=l k 

where rj represent the position vectors of the N constituent electrons (or holes); aj,ak 
denote the creation and annihilation operators for LO phonons with wave vector k. The 
/k were determined variationally [73]. It may be emphasized that the expression l|21|) 
constitutes - especially in its implementation - a nontrivial extension of the one-particle 
approximation to a many-body system. An advantage of the many-polaron canonical 
transformations introduced in Ref. (73] for the calculation of the ground state energy of 
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a polaron gas is that, to order a, the many-body effects are entirely contained in the 
static structure factor of the electron (or hole) system, which appears in the analytical 
expression for the energy. The minimum of the total ground-state energy per particle as 
a function of the electron density for a polaron gas is shown [73] to lie at lower value of 
the density than that for the electron gas. 

In Ref. |25j. the optical absorption coefficient of an interacting many-polaron gas 
has been derived starting from the many-polaron canonical transformations and the 
variational many-polaron wave function introduced in Ref. |73| . The real part of the 
optical conductivity of the many-polaron system to order a is obtained in terms of the 
dynamical structure factor in Ref. [35] in the form which has been discussed in the 
Subsection 12' II The experimental peaks in the mid-infrared optical absorption spectra 
of fFig. HUf) and manganites fFig.lll[) have been interpreted within this theory, using 
0. The experimental optical absorption spectrum (up to 3000 cm -1 ) of Nd 2 Cu02-5 




Fig. 10. - The infrared absorption of Nd2Cu02-<s (5 < 0.004) is shown as a function of frequency, 
up to 3000 cm -1 . The experimental results of Calvani and co-workers |82| is represented by 
the thin black curve and by the shaded area. The so-called 'd-band' rises in intensity around 
600 cm -1 and increases in intensity up to a maximum around 1000 cm -1 . The dotted curve 
shows the single polaron result. The full black curve represents the theoretical results obtained 
in the present work for the interacting many-polaron gas with no = 1.5 x 10 17 cm -3 , a — 2.1 
and mb = 0.5 m e . (Recently, a numerical inaccuracy was identified in the fitting procedure, 
concerning the strength of the plasmon branch, so that the parameter values given above may 
not be the optimal parameter values, see Ref. 83 .) (From Ref. |25|.) 



(5 < 0.004), obtained by P. Calvani and co-workers |82|. is shown in Fig. HOIfshaded area) 
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Fig. 11. - Comparison of the measured mid-infrared optical conductivity in La 2 /3Sr 1/ / 3 Mn03 
at T = 6 K to that given by several model calculations for m = 3m c , a of the order of 1 and 
no = 6 x 10 21 cm -3 . The one-polaron approximations [the weak-coupling approach by V. L. 
Gurevich, I. G. Lang, and Yu. A. Firsov [24) (GLF model) and the phenomenological approach 
by D. Emin [841 (E model)] lead to narrower polaron peaks than a peak with maximum at 
uj ~ 900 cm -1 given by the many-polaron treatment by J. Tempere and J. T. Devreese (TD 
model) of Ref. [21]. (After Ch. Hartinger et at, Ref. pi].) 



together with the theoretical curve obtained by the method of Ref. |25] (full, bold curve) 
and, for reference, with the one-polaron optical absorption result (dotted curve). At 
lower frequencies (600-1000 cm -1 ) a marked difference between the single polaron optical 
absorption and the many-polaron result is manifest. The experimental rf-band can be 
clearly identified, rising in intensity at about 600 cm^ 1 , peaking around 1000 cm -1 , and 
then decreasing in intensity above that frequency. At a density of no — 1.5 x 10 17 cm~ 3 , 
we found compelling agreement between our theoretical predictions and the experimental 
curve. The authors of Ref. [SS] interpret their data on La2/3Sri/3Mn)3 with Eq. Q 25 , 
adapted for an on-site Hubbard interaction. 

The optical conductivity of a many-polaron gas was further investigated, by Iadonisi 
and associates, e. g. in Ref. 86 within the memory-function approach |.'i5| by calculating 
the correction to the dielectric function of the electron gas, due to the electron-phonon 
interaction with variational parameters of a single-polaron Feynman model. A suppres- 
sion of the optical absorption from the one-polaron optical absorption of Ref. [35] with 
increasing density shown in Fig. ^]is expected because of the screening of the Frohlich 
interaction with increasing polaron density. 

Recently, dynamical exchange effects were found to lead to noticeable corrections to 
the random-phase approximation results for the ground state energy, the effective mass, 
and the optical conductivity of a two-dimensional many-polaron gas in Refs. [871 183j . 
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Fig. 12. - Optical conductivity of a polaron gas at T = as a function of the frequency as 
calculated in Ref. )8fci| (CDI) for different values of the electron density: no = 1.4 x 10 (solid 
curve), no = 1.4 x 10~ 4 (dashed curve), no = 1.4 x 10 -3 (dotted curve), and no = 1.4 x 10~ 2 
(dotted-dashed curve). The electron density is measured per Rp, where R p is the Frohlich 
polaron radius. The value of £o/£oo is 3.4. The solid curve practically coincides with the known 
optical conductivity of a single polaron 35 (DDG). (After Ref. |86|.1 



4"2. Many-Polaron Cyclotron Resonance in Quantum Wells. - In GaAs/AlAs quan- 
tum wells with high electron density, anticrossing of the cyclotron-resonance spectra has 
been observed near the GaAs transverse optical (TO) phonon frequency ujti rather than 
near the GaAs LO-phonon frequency &li (88| . This anticrossing near u>ti was explained 
in the framework of the polaron theory in Refs. (Fig. ll3|) . This effect is in con- 

trast with the cyclotron resonance of a low-density polaron gas in a quantum well, where 
anticrossing occurs near the LO-phonon frequency. The appearance of the anticrossing 
frequency close to u>ti instead of wja is due to the screening of the electron-phonon inter- 
action by the plasma oscillations. The theory takes into account the magnetoplasmon- 
phonon mixing and the electron-phonon interaction with all the phonon modes specific 
for the quantum well. A discussion on alternative explanations of the data of Ref. |HH1 
is still going on [9U1 13T1 192] . 

4'3. Interacting Polarons in a Quantum Dot. - A system of N electrons (or holes) with 
mutual Coulomb repulsion and interacting with the bulk phonons is analysed in Refs. 
I64| . It has been shown |M] that the path-integral approach to the many-body problem 
for a fixed number of identical particles can be formulated as a Feynman-Kac functional 
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Fig. 13. - Density map of the magnetoabsorption spectrum for a 10-nm GaAs/AlAs quantum 
well as calculated in Ref. |89| . Peak positions of the experimental spectra are taken from Fig. 3 
of Ref. |88|. Crosses indicate the CR peak positions for B < 17 T. Above 17 T, the CR line splits 
into two components represented by triangles and by filled dots. This splitting occurs due to the 
fact, that the electron Landau levels in a non-parabolic conduction band are non-equidistant. 
The higher- and lower-energy components correspond to the transitions between the Landau 
levels (0 — * 1) and (1 — > 2), respectively. A single CR line shown by squares is observed above 
23 T. Dashed lines show LO- and TO-phonon energies in GaAs. (After Ref. |89|.^ 

on a state space for N indistinguishable particles. The resulting variational inequality for 
identical particles formally has the same structure as the Jensen- Feynman variational 
principle [TOj- A model system consisting of N electrons and Nf fictitious particles 
in a harmonic confinement potential with elastic interparticle interactions is studied 
in Refs. |931 I64j . The Lagrangian of this model system is symmetric with respect 
to electron permutations. The parameters of the model system are found from the 
variational principle 95; and used for the calculation of the ground-state energy and the 
optical conductivity of the iV-polaron system in a quantum dot 1^ . 

The polaron contribution to the ground-state energy of the iV-polaron system in a 
quantum dot per particle calculated within our variational path-integral treatment for 
different N |97) is compared with the respective quantity for a polaron gas in bulk, 
found in Ref. |99j within an intermediate-coupling approach and in Ref. jlUOj . using a 
variational approach developed first in Ref. |73| . Our all-coupling variational method 
|97j provides somewhat lower values for the polaron contribution than those obtained in 
Refs. 123 [100j and the trends predicted in those references are confirmed. 
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The "symmetric" polaron model (Devreese 1964) 
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Fig. 14. - Upper panel: The phonon distribution functions Wat in the "symmetric" polaron 
model for various values of the effective coupling constant A at x = 1,P = (from Fig. 
23). Lower panel: Distribution of multiphonon states in the polaron cloud within DQMC 
method for various values of a (from |23| . Fig. 7). 
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Appendix A. 

On the Contributions of the A-phonon States to the Polaron Ground State 

In Ref. - the present author introduced and analysed an exactly solvable ("sym- 
metric") lD-polaron model. The further study of this model was performed in Refs. 
|32l 1101) . This model can be used as a tool for the evaluation of the accuracy of the 
standard polaron theories [TU51 11051 HTM] . 

The model consists of an electron interacting with two oscillators possessing opposite 
wave vectors: k and — k. The parity operator, which changes ak and a_k (and also a k 
and ilk)) commutes with the Hamiltonian of the system. Hence, the polaron states are 
classified into even and odd states with eigenvalues of the parity operator +1 and — 1, 
respectively. For the lowest even and odd states, the phonon distribution functions Wn 
are plotted in Fig. 1141 upper panel, for some values of the effective coupling constant 



A of this "symmetric" model. The value of the parameter x = y (hk) /mjfewLO for 
these graphs was taken to be 1, while the total polaron momentum P = 0. In the 
weak-coupling case (A « 0.6) Wn is a decaying function of N. When increasing A, Wn 
acquires a maximum, e. g. at N = 8 for the lowest even state at A = 5.0625. The phonon 
distribution function Wn has the same character for the lowest even and the lowest odd 
states at all values of the number of the virtual phonons in the ground state (as distinct 
from the higher states). This led to the conclusion that the lowest odd state is an internal 
excited state of the polaron. 

In Ref. the structure of the Frohlich polaron cloud was investigated using the 

DQMC method. Contributions of 7V-phonon states to the polaron ground state were 
calculated as a function of N for a few values of the coupling constant a, see Fig. 1141 
lower panel. The evolution from the weak-coupling case (a — 1) into the strong-coupling 
regime (a = 17) was studied. Comparison of the lower panel to the upper panel in Fig. 
EJshows that the evolution of the shape and the scale of the distribution of the iV-phonon 
states with increasing a as derived for a Frohlich polaron within DQMC method [S3] is 
in notable agreement with the results obtained within the -'symmetric" ID polaron model 
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